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An eikonal expansion is developed in order to provide systematic corrections to the eikonal ap- 
proximation through order 1/fe 2 , where k is the wave number. The expansion is applied to wave 
functions for the Klein-Gordon equation and for the Dirac equation with a Coulomb potential. Con- 
vergence is rapid at energies above about 250 MeV. Analytical results for the eikonal wave functions 
are obtained for a simple analytical form of the Coulomb potential of a nucleus. They are used to 
investigate distorted-wave matrix elements for quasi-elastic electron scattering from a nucleus. Fo- 
cusing factors are shown to arise from the corrections to the eikonal approximation. A precise form 
of the effective-momentum approximation is developed by use of a momentum shift that depends 
on the electron's energy loss. 

PACS numbers: 



I. INTRODUCTION 

Quasielastic scattering of electrons by nuclei can pro- 
vide important information about the response of a nu- 
cleus to a weakly interacting probe. Experiments have 
been performed at the MIT Bates Laboratory [1-9], at 
the Saclay Laboratory [10-15] and at SLAC [16-18] in 
order to explore this reaction. A recent review of quasi- 
elastic scattering provides a guide to the experimental 
and theoretical results [19]. One important theoretical is- 
sue concerns the corrections that arise from the Coulomb 
potential of the nucleus. Another issue is the accuracy 
with which longitudinal (L) and transverse (T) response 
functions can be determined when Coulomb corrections 
are present. 

In order to include the Coulomb distortion of electron 
waves, it is necessary to solve the Dirac equation for scat- 
tering of an electron of mass m and energy E from a 
nucleus, i.e., 



E - a ■ p - V(r) - (3m)>(p = 



(1) 



where V(r) is the Coulomb potential due to the nu- 
clear charge distribution. Exact solutions for the Dirac- 
Coulomb wave functions may be obtained as a sum over 
partial waves. [20-22] However, the sum involves deli- 
cate numerical computations. As the electron energy 
increases, the partial- wave expansions converge more 
slowly in spite of the fact that the Coulomb corrections 
become smaller. A simpler approach is to use the eikonal 
approximation for the distortions because the approxima- 
tion gets better as the energy increases and it can provide 
insight into the parameters that control the Coulomb cor- 
rections. 

For waves moving along the z-direction, the eikonal 
approximation provides an approximate solution, 



where k = ^/EP- 



is the momentum and \ is a sim- 



e ikz e i X 



(2) 



pie function of the potential. It is important to include 
focusing factors such as those derived by Yennie, Boos 
and Ravenhall [23-25] based upon the WKB approxi- 
mation to the distorted Coulomb waves. Some partic- 
ulary transparent results have been obtained using the 
eikonal approximation to derive an effective-momentum 
approximation (ema) [26, 27] that produces results very 
similar to plane-wave results. However, with the combi- 
nation of the eikonal approximation, WKB focusing fac- 
tors and the effective-momentum approximation, the ap- 
proach lacks a systematic basis and accuracy is uncertain. 
Approximations to the partial-wave analysis, such as the 
analysis of Ref [28] , that attempt to improve the effective- 
momentum approximation also have uncertain accuracy. 
Significant disagreements in the determination of nuclear 
response functions from experimental data [19, 29] have 
arisen at least in part owing to the use of different theo- 
retical methods to remove the Coulomb corrections. 

In order to address the issue of Coulomb correc- 
tions, we develop corrections to the eikonal approxima- 
tion based upon a systematic expansion about the high- 
energy limit. This eikonal expansion is shown to be 
rapidly convergent for typical energies and targets used in 
quasi-elastic scattering. For a simple analytical Coulomb 
potential, analytical forms are developed for the eikonal 
corrections through second-order in 1/k, where k is the 
electron's wave number. 

The eikonal approximation has a long history [30]. 
Glauber used the eikonal approximation to develop a sim- 
ple and insightful form of multiple scattering theory that 
is called Glauber theory [31]. Czyz and Gottfreid [32] 
used the eikonal approximation to analyze electron scat- 
tering but that work did not include focusing factors. 
Work by Giusti et al. also is based on the eikonal ap- 
proximation [33, 34] and some recents works have com- 
bined the eikonal approximation with semi-classical fo- 
cusing factors in order to assess Coulomb corrections in 
quasi-elastic scattering. [35, 36]. We show that the fo- 
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cusing factors arise from systematic corrections to the 
cikonal approximation. 



Corrections to the eikonal approximation also have a 
long history. Work by Saxon and Schiff [37] showed how 
to correct the approximation to leading order in 1/fc. A 
systematic expansion for the scattering t-matrix was de- 
veloped by Sugar and Blankcnbecler [38] . Systematic cor- 
rections to the Glauber approximation were developed by 
Wallace [39] and extended to the Dirac scattering ampli- 
tude in [40]. However, a systematic expansion for wave 
functions has not been developed prior to this work. 



In Section II, we develop the eikonal expansion for 
Klein-Gordon wave functions because that is the sim- 
plest and most transparent case. In Section III, we de- 
velop the eikonal expansion for the Dirac wave function. 
Section III A focuses on u(r), which is a Pauli spinor 
containing the two upper components of the Dirac wave 
function. Because there is a spin-orbit interaction, addi- 
tional spin-dependent terms arise in the eikonal expan- 
sion. Because the analysis is technically more compli- 
cated, details are given in an appendix. In Section IIIB 
we show that the Pauli spinor £(r) that contains the two 
lower components of the Dirac wave function takes a very 
simple form in the limit of vanishing electron mass, i.e., 
£(r) = 2Xu(r), where A = ±| is the helicity. This fol- 
lows because of the structure of the Dirac equation in 
the limit of vanishing electron mass. Conservation of the 
electron's helicity holds to high accuracy at electron ener- 
gies of interest. Section III C discusses convergence of the 
eikonal expansion and analytical results for the eikonal 
phases, which are given in an appendix. Section III D dis- 
cusses the different focusing factors that arise for Klein- 
Gordon and Dirac waves. In Section IV, we consider 
helicity matrix elements of the electron current based 
on the Dirac distorted waves. We show that the spin- 
orbit distortions arising from the Coulomb interaction 
can alter the longitudinal and transverse helicity matrix 
elements in somewhat different ways, which can affect 
the L/T separation. Specializing to the longitudinal re- 
sponse, section V discusses quasi-elastic scattering by use 
of a simple model of the nuclear response. Section VI re- 
visits the effective-momentum approximation (ema) and 
Section VII presents results of calculations of the longi- 
tudinal response function. The analytical eikonal wave 
functions with systematic corrections included are used 
to describe the Coulomb corrections. We show that the 
effective-momentum approximation can be made precise 
by use of a calculated momentum shift that depends on 
the electron's energy loss, 10. A summary and some con- 
clusions are presented in Section VIII. 



II. EIKONAL EXPANSION FOR 
KLEIN-GORDON WAVES 

For the Klein-Gordon equation with a Coulomb poten- 
tial V(r), one has 

{[E-V(r)] 2 -p 2 -m 2 )^(r)=0. (3) 

The potential may be a point-like Coulomb potential or a 
potential that is derived from a finite charge distribution, 
such as that of a nucleus. Writing the wave function in 
the form of a plane wave propagat ing in the z-direction 
with wave number k = \JE 2 — m 2 and a complex phase 
shift x, 

V>(r) = e 4fc V x(r) , (4) 

leads to the following equation for the phase shift, 

[E - V{r)] 2 - [kz + Vx - iV] • [kz + W X ] - m 2 = 0. (5) 

Using the fact the E 2 — k 2 — m 2 = and dividing by 2k 
leads to a differential equation for x, 

d x _ v v 2 (Vx) 2 ,v 2 x (6) 

dz v 2k 2k 1 2k ' U 

where v = k/E m 1 when E » m. 

For outgoing-wave boundary conditions and a poten- 
tial that decreases faster than 1/r, the wave function 
must be a pure plane wave as z — > — oo and a phase- 
shifted plane wave as z — > oo. For a Coulomb potential, 
the same boundary condition is used with the under- 
standing that the potential is cut off at a large distance 
r > A. For either case, the outgoing-wave eikonal phase 
vanishes as z — > — oo and the eikonal phase is found by 
integrating Eq. (6) as follows, 

X i+) {v) = -- f dz'V{r') 

v J —oo 

-Yk /^'(W +) (r')) 2 -^) 

+ Yk fj^ l2 x (+ \*')- (7) 

A superscript (+) denotes the outgoing-wave boundary 
condition. 

The eikonal expansion is an iterative solution of Eq. (7) 
about the limit k — > oo. It initially takes the form 

X i+) =X ( o +) +x[ +) +X ( 2 +) +---, (8) 

where the "barred" phase shifts are complex and their 
subscripts denote the order of iteration. The lowest order 
term is appropriate to the limit k — > oo and is obtained 
from the first term on the right side of Eq. (7), 

X ( +) (r) = -i f dz'V(r'), (9) 

u J —oo 
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where r' = V z' 2 + b 2 and b is the impact parameter. 

Higher order terms in the expansion are evaluated as 
follows. In first order, one evaluates the left side of 

Eq. (7) using + Xi + ^ an d the right side using Xo^- 
After cancelling terms that are equal because of Eq. (9), 
the first-order term xi + ' ) is obtained as 



s/_> s '( <v '*» +>(r ' ,)2 - l/2 ) 

+ lk fj h 'V\V\*')- (io) 

In second order, one evaluates the left hand side of Eq. (7) 
using Xg + ^ + x! + ' + %2 + ' an d the right hand side using 
Xo + ^ +Xi + ' ) • After cancelling terms that are equal because 
of Eqs. (9) and (10), the second-order term is obtained 
as 



x' 2 + > 



dz' 



+ (V'xi +) (r')) 2 



+ 



2V'^ +) (r')-V'x< +) (r') 

— / rf,A7' 2 ,-(+)^ 
2fc 



^'V'xS + '(r'). 



(11) 



Higher-order terms are not considered in this work. 

At the final stage of analysis, it is convenient to sepa- 
rate the real and imaginary parts of the eikonal phase as 
follows, 



X {+) = X (+) + 



(12) 



and to expand each of these "unbarred" phases in sys- 
tematic powers of 1/fc, 



X 



(+) 
(+) 



x[ +) +xi +) 



- U), 



(+) 



(13) 



Here the subscripts denote the power of 1/fc that is in- 
volved. Through second order, the systematically or- 
dered phases for the Klein-Gordon waves are, 



X^ +) (r) = 



Xi +) W = 
X^ +) (r) = 
+ 



dz'V(r') 

k / 2 oo ^'( (v ' x ° +)(r ' ))2 - y2(r ' ) ) 



v'M + V) 



2V'x^ +) (r')-V'xr(r') 



4 +) (r) = ± f Jz'V' 2 X { +\v') 



1 

2fc 



dz' 



2V'x^ +) (r')-V'^(r') 



(14) 



For the Klein-Gordon wave function with outgoing- 
wave boundary conditions, this gives 



^ (+) (r) =e tkz e ix e 



(+) _,.,(+) 



(15) 



and one may work at various orders in the eikonal expan- 
sion by truncating the expansions of Eq. (14). 

Incoming wave boundary conditions are appropriate 
for final-state wave functions in matrix elements. In that 
case the wave function must be a pure plane wave as 
z — > oo and a phase-shifted plane wave as z — > — oo. It 
is written as 

^-\r)=e ikz e- i * ( ~\ (16) 

which leads to 

dx (~) _ V V 2 (Vx ( -)) 2 , .v 2 x<-> 



dz 



v 2k 



+ 



2k 



2k 



(17) 



Integration with incoming-wave boundary conditions 
produces 



X (_) (r) 



dz'V(r') 

/■OO 

-j dz'V 2 X (-\r'). I 



i 

2k 



As before, the eikonal expansion is an iterative solution 
of Eq. (18) that produces an expansion about the limit 
fc — > oo. It initially takes the form 



X (_) =Xo ] +x[ ] +X { 2 ] + ■■ 
with the lowest order term being 



X^ 



(r) = - l -j°° dz'V{r'). 



(19) 



(20) 



Higher order terms in the expansion are evaluated in 
the same manner as discussed above. They are 



Xi" } (r) = 



X ( 2~\r) 



~k\ z dz 'V' 2 X { -\r'), 
~2kJ z 



2V'x^ ) (r')-V'xi" ) (r') 



dz'V 2 xf\v'). 

(21) 

Separating the real and imaginary parts of the eikonal 
phases as follows, 



X ( ) = X ( ) - «w ( \ 



(22) 
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and expanding each of these "unbarred" phases in sys- 
tematic powers of 1/fc, we have 



X ( } = Xo } + Xi ] + x[ ] + ■ ■ 



,(-) 



,,(-)4-,,(-) 4- 



(23) 



Such systematically ordered phases for the Klein-Gordon 
waves for incoming wave boundary conditions may be 
obtained from Eq. (14) by replacing all integrations /_ 

by and all superscripts (+) by (— ). Because of the 
symmetry of the potential with respect to inversion of z, 
the following relations hold for each order n, 



xi-Hz,b) = x ( +H-z,b), 

^-\z,b)=J+\-z,b). 



(24) 



For the Klein-Gordon wave function with incoming wave 
boundary conditions, this gives 

</,(-) ( r ) = e ikz e- ix( ~ ) e-^~\ (25) 

and one may work at various orders by truncating the 
expansions for the eikonal phases. 

The imaginary part of the eikonal phase produces the 
"focusing factor" , 



f KG(±) {r) = e - 



u,<±>(r) 



(26) 



In order to provide some insight into the focusing factor, 
it is useful to consider a simple potential for which the 
eikonal phases may be determined analytically. That is 
done in Appendix A and we find that at r = 0, 



(± )(Q) _ V(0) „ V(0) 



(27) 



where V(Q) is the potential at the origin. Thus, the fo- 
cusing factor is approximately 



f KG W(0)tvl-V(0)/(2E). 



(28) 



A matrix element for emission of a photon with mo- 
mentum q and energy u> = Ei — Ef involves initial and 
final momenta and kj. It takes the form 

M " = J dMi- ] *{v)f KG e-^^\v), 

= f Sre^-^e^ff G '-\v)fr (+ \r)f K m 



where Q = k^ — k f and \ — X/ ' + Xi + ^ ■ It should be 

noted that x\ an d ^\ are obtained from Eqs. (14) 
with the z-axis parallel to initial momentum k^, while 
' and ojf ' are obtained from the same equations 

using the outgoing- wave condition (/°°) with the z- 
direction parallel to final momentum k f. 



The conserved current for the Klein-Gordon equation 
also contains Coulomb distortions, i.e., 



3kg - [Ei + E f - 2V(r)] I^E~E S , 
juG = [ki + k f + Vx< +) (r) - Vx^(r)] / 1 ^E~E h 

(30) 

where phase-space factors for initial and final states have 
been included. Current conservation holds because 

^3 kg - q • 3K G )e-^ r cx e-^G-^Ei) - 

G-\E f )e-^\ (31) 

where G^ 1 (E) = [E — V(r)] 2 — p 2 - m 2 annihilates the 
wave function at energy E. Relative to the plane-wave 
current, the Klein-Gordon current j® KG contributes to the 
focusing effect a factor equal to (1 — V(r)/E) where E = 
\{Ei + Ej) is the average energy. The initial-state wave 

function provides a factor given by Eq. (26), i.e., / i ° < ' + ' ) ~ 
(1— V(0)/ (2E)) and the final state wave function provides 
a similar factor. The combination of these factors at 
r = is approximately equal to (1 — V(0)/E) 2 . That 
agrees with the result of Yennie et al. [23] for the Dirac 
wave function. 

In the case of Dirac-Coulomb waves, the conserved 
electron current is the Dirac matrix 7 M , which it is not 
modified by the presence of a Coulomb potential. It does 
not contain any focusing effects. The Dirac focusing fac- 
tors arise solely from the wave functions and are given 
in Eq. (51), which yields / D (±) w (1 - V(0)/E) for each 
wave function. Thus, there is a different focusing fac- 
tor than for a Klein-Gordon wave function. However, 
the combination of currents and wave functions produces 
similar overall focusing effects for the Klein-Gordon and 
Dirac current matrix elements. We note that the analy- 
sis of Refs. [35, 36] used the Dirac focusing factors with 
the Klein-Gordon current, which produces one too many 
factors of 1 - V(0)/E. 



III. EIKONAL EXPANSION FOR 
DIRAC-COULOMB WAVES 

For the Dirac equation, the eikonal expansion is carried 
out in two stages. First we consider the Pauli spinor u(r) 
that contains the two upper components of the Dirac 
wave function, i.e., 



V(r) 



u(r) 



(32) 



It follows from the Dirac equation that the Pauli spinor 
£(r) that contains the two lower components may be de- 
termined in a second stage from 



£(v) 



where E-i — E + m. 



1 



En- V 



a ■ p u(r), 



(33) 
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A. Upper component spinor 



Eliminating the lower component spinor from the 
Dirac equation leads to the following equation for the 
upper-component spinor 



E 1 -V-a-p 



E 2 -V 



a ■ p Wr), 



(34) 



where E\ = E — m. For electron scattering it is generally 
the case that E » m and thus Ex m E 2 ~ E. 

For outgoing-wave boundary conditions, the Pauli 
spinor u(r) is written in terms of a complex eikonal phase 
X^ + H r ) an d a complex spin-dependent phase 7^ + )(r) as 
follows 



(r) = 1 



V 
E 2 



1/2 



ikz ix (+> i(7el {+) 



(35) 



The wave propagates in the ^-direction and an impact 
vector b is defined as the part of r that is perpendicular 
to the z-direction, i.e., b = z x (r x z). Three orthogonal 
unit vectors are : z, b = b/|b| and e — by, z. The spin 
matrix in the eikonal phase is a e — a ■ e. The factor 
(1 — V/E2) 1 / 2 is introduced in order to sum up terms 
that otherwise arise in higher orders. 

The eikonal expansion is developed in Appendix B. 
The result is that the eikonal phases are expanded in a 
systematic fashion in powers of 1/k as follows. 



x (+) = 



x i +) +x[ +) +x 2 +) + 



(+) = 
(+) - 



,(+) 



+ LO. 



(+) 



7 



+ 
+ ■ 



4 +) + ■ 



(36) 



for the Dirac equation are as follows 



x[ + \r 

x[ + \ 



7^ +) (r 



dz'V(r'), 

= -Yk f Jz'^'x^ {*')?- V\r>)), 
(2V'x^ +) (r')-V'xl +) (r') 



1 

2k 



dz' 



+ v'M + V) 



AkE 2 



= Yk /^'[v'M+V) 

- 2V' X ( +) (r')-V'^ +) (r')" 
2k L J db 



1 

2k 



dz'\2Vx^\v')-V 1 [ + \v') 



(+)' 



1 dVdx'p 
E dz db 



#'v-hL"{ v '-*y +>V) - (37) 

Results for Xo + \ Xi + \ an d t»4 + ^ are the same as for 
the Klein-Gordon case. The last term in the result for 
X2 + ^ is not present in the Klein-Gordon case. 

Spin-dependent phases can be shown to be simply re- 
lated to the spin-independent ones as follows, 



1 



d 



E 2 -V db 



This connection holds to order 1/k 2 . 

The upper-component spinor of the Dirac wave func- 
tion for hclicity A and outgoing-wave boundary condi- 
tions is given by 



,(+) 



V \V2 



(38) 



where £a is a helicity eigenstate. One may work at var- 
ious orders of the eikonal expansion by truncating the 
expansions of Eq. (36). 

The upper-component spinor for helicity A and 
incoming-wave boundary conditions is given by 



e~" 1 ) ^ (39) 



where 



where the subscript of each term denotes the power of 
1/k that is involved. The systematically ordered phases 



7<") - iS^ 
7^+7^ 



+ ...), 



(40) 
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with phases obtained from Eqs. (36) and (37) by re- 
placing the superscripts (+) by (— ) and the integra- 
tion ranges J* by Alternatively, the symmetry 
of Eq. (24) may be used. 



B. Lower-component spinors 

Equation (33) specifies the lower-component spinor in 
terms of the upper-component spinor. The connection is 
very simple in the limit m — > when helicity eigenstates 
are used. We find 



h(r) = 2Xu x (r) 



C. Convergence of eikonal expansion 

Convergence of the eikonal expansion is discussed for 
scattering of a 500 MeV electron by 208 Pb. Given that 
the electron mass is m = .511MeV, it follows that k s=s E 
and »wl, both within a part per million. The Coulomb 
potential is approximately 1^(0) = 25 MeV at the cen- 
ter of the nucleus. The leading order eikonal phase 
is of order 2RV(0), where R is the mean radius and 
the factor 2R provides an estimate of the integral over 
z. If the charge distribution is approximated as con- 
stant within a sphere of radius R, then V(0) = ^Za/R, 

^42) so we expect that ~ iZa. For 208 Pb, this gives 



where A = ±|. This holds generally in the m = limit 
as the following analysis shows. 

For m = 0, we have E\ = E 2 — E. If we define 

/ \ 1/2 

u\= (I — V/E\ u\(r), then Eq. (34) may be written 

as 

u A (r) = h 2 u x , (42) 
where for E > \ V\, h is the hermitian operator 
1 I 



h = 



VE-V 



G ■ p 



y/E-V' 



(43) 



Thus, u\ is an eigenfunction of the hermitian operator h 2 
with eigenvalue +1. It must then also be an eigenfunction 
of h with eigenvalues ±1, that is, 

hu x = ±u x . (44) 

However, because of the Dirac equation it follows that 

hux = ix, (45) 

1/2 



where ix 
that 



(i-v/e) 



ix ■ These two equations require 



lx = 2Xux 



(46) 



where A = ±|, which proves Eq. (41). Clearly, h is the 
helicity operator for the m = Dirac equation. Con- 
servation of helicity is a well-known result for the limit 
m — > when the interaction is a vector current, or 
any single component of a vector interaction as for the 
Coulomb interaction. In the limit as V — > 0, h becomes 
the usual helicity operator for a plane-wave state. 

For a nonzero electron mass, a similar but approximate 
analysis may be performed to show that 




E 2 - (V) . 



(47) 



where (V) is the average potential. The correction term 
is about one part per thousand for a 500 MeV electron. 
It will be neglected in this work. 



Xq + ^ ~ 2. The eikonal expansion introduces corrections 



that involve the nondimensional ratio V/E w .05, so we 

expect xi +) ~ -0025xo +) « -005. The eikonal expan- 
sion is not convergent but is asymptotic, meaning that 
the error should be bounded by the first neglected term. 
When terms up to second order are kept, the error is 
of order Xa +) ~ .000125xo +) « .00025 in the example 
discussed here. There is a second expansion parameter 
involved in terms like u>\ or u) 2 , namely 1/ka, where a 
is a length parameter that characterizes derivatives of 
the potential, i.e., |VV(r)| w l/a\V(r)\. Provided that 
the potential is sufficiently smooth, these corrections also 
are small. Thus, the expansion can produce accurate 
wave functions for electron scattering when the energy is 
sufficiently high and the potential is sufficiently smooth. 
One must be aware that the eikonal wave functions omit 
the wave bending effects that arise from exact solutions 
and are less accurate far from the nucleus. However, the 
eikonal wave functions can provide an accurate descrip- 
tion of Coulomb distortion effects in the region where the 
nuclear charge density is nonzero. 

It is possible to reduce the terms in the eikonal expan- 
sion to analytical forms for the potential 



V(r) = 



aZ 



Vr 2 + R 2 



(48) 



where a = e 2 /(hc), Z is the nuclear charge and R is a 
length parameter. This Coulombic potential corresponds 
to a charge density 



p(r) = 



3Ze 



I? 2 



4tt (r 2 + i? 2 ) 5 /2 ' 



(49) 



The analytical results are given in Appendix A. 

Figure 1 shows the eikonal phases for a charge Z = 100, 
electron energy E = 200 MeV and radius parameter R = 
2 fermi. These parameters are chosen in order to make 
the corrections visible. The corrections are much smaller 
for a 500 MeV electron and smaller nuclear charge. 



D. Focusing factors 

As noted above, the focusing factors are important 
when Coulomb distorted waves are used. Ignoring for 



7 




b(F) 



FIG. 1: Eikonal phases at z=0 versus impact parameter: solid 
line shows xo, dashed line shows xi, dotted line shows X2, 
rectangles show uii and ovals show ui2- A constant has been 
added to xo such that it vanishes at b — 3.5R. Phases are 
shown for Z=100, E=200 MeV and R=2 fermi. 



as may be seen by comparing Eqs. (26) and (51). The 

e~ w i +) f« 1 — V(0)/(2Ei) factor is the same to leading 
order for a Klein-Gordon and a Dirac wave function. 
Combining this part of the eikonal correction with the 
(1 - V/E 2i ) 1/2 yields a focusing factor w 1 - V/E l 

in the Dirac wave function, thus reproducing at r = 
the expected factor 1 — V(0)/Ei that has been derived 
by Yennie, Boos and Ravenhall [23] based on a WKB 
analysis of the Dirac-Coulomb wave function. A simi- 
lar result holds for the final-state focusing factor, \ 
which is approximately 1—V/Ej. Thus, the overall focus- 
ing effect in the matrix clement is approximately equal to 
(1-V(0)/E f )(l-V(0)/Ei), which is similar to the over- 
all effect in the Klein-Gordan matrix element, Eq. (29). 
The reason that focusing factors did not appear in early 
analyses based on the eikonal approximation is because 
they arise from corrections to the eikonal approximation. 

In passing, we note that the Glauber approximation 
is obtained when the eikonal phases for initial and fi- 
nal states are evaluated using for each a z-axis paral- 
lel to the average momentum, + kf), and only the 

leading-order phases, Xo + '' an d Xo \ are retained. This 
approximation omits the focusing factors. 



the moment the lower component spinors, we consider a 
current matrix element for emission of a photon of energy 
lu = Ei — Ef and momentum q using upper-component 
spinors corresponding to initial momentum k^, initial hc- 
licity Xi, final momentum kf and final hclicity A/. The 
momentum transfer is Q = kj — kj and the current ma- 
trix clement is 



M" = / d 3 ru^~ ) *(r)7"e- iq - p t*^(r) 



(+), 



J d 3 rCi / e l(Q - q) - r e^e^/^ ) */f- ) 7^/ J D(+) x 



v< + ) 



,(-)/ 



(50) 



where x — X/ ; ( r ) + Xi + ^ ( r ) an d the Dirac focusing fac- 
tors are defined by 



/f (+) (r) = 1- 



\ 1/2 

E 2 J 6 ' 



f? { -^= (l-^) 1/a e^-\ (51) 



Note that \\ 



(+) ,,(+) 



, and 7 4 (+) = 7 2 (+) + iS^' are ob 



(+) 



tained from Eqs. (36) to (37) with the z-axis parallel 

to initial momentum k^, while X/~\ w i \ an d l\ ' = 

7^ — iS\ ' , are obtained from the same equations using 
the outgoing-wave condition (J ) with the z-direction 
parallel to final momentum ky. 

Focusing factors differ from those appropriate to a 
Klein-Gordon wave function by the factor (1 — VjE^i) 1 ^ 2 



IV. ELECTRON CURRENT MATRIX 
ELEMENTS 

Electron scattering involves the current matrix element 

ft = J d 3 r^-{ f (r)re-^i + X (r), (52) 

where q is the three-momentum of a photon emitted at 
point r, uj = Ei — Ef is the energy of the photon and 



T (+) , s 1 f u { +) ( r )\ 



(53) 



is a Dirac-Coulomb wave with a normalization factor 
1 / \pl included in order that it reduces as V — > to the 
m = plane- wave spinor, 



1 

72 



1 

a ■ k, 



(54) 



Using lower-component spinors from Eq. (41) and the 
conventions of Bjorken and Drcll [41] for the 7 M matrices, 
one readily finds that the matrix element j® involves the 
overall factor 1 + 4A/Aj = 28\.\ i and the matrix element 
j e involves the overall factor 2A/ + 2A, = (2\ i )25\ f \ i . 
Thus, helicity is conserved as it must be in the m = 
limit. The electron current matrix elements are further 
reduced by use of the convention of Kubis [42] for the 
helicity matrix elements between initial and final states, 
leading to 



= i(Q-q)-r ix f D (~h 



D(+), 



(55) 
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where h%(r) is a four- vector of helicity matrix elements, 



,-,(-)* 



,-,(-) 



{h° e ,h e } = C{ f (0e)e^f '~{l,a}e^ 



-(-)* 

COSJj: 



cosj^ + ia ei siwyl^ 



r sinjj: 



(-)- 



(56) 



Here £ x (0 e ) = ^^ e J<J y e/ 2 i s the helicity eigenstate for 
the outgoing electron and e is the scattering angle of the 
electron. We find that 



(57) 



Angle <t>i is the polar angle of initial state impact pa- 
rameter, i.e., bi = cosSii + sincfuy and Sf is the polar 
angle of the final state impact parameter, where initial 
and final momenta are in the xz-plane: = kiZ and 
kf = kf[cos0 e z + sinO e x\. Moreover the required helicity 
matrix elements are 

£\ f {8e)Z.\i = $\ f \iCos-6 e + (Xf - \i)sin-9 e 

£\ f (0e) <? = (Xf + Xi)(e 2 XiSin^e e + i z cos^6 e ^J 

+ |A/ - Ai| (i 2Xi cos^9 e - e z sin^e e ^j (58) 



where e 2X = e x + (2X)ie y . 

Combining these parts produces the required helic- 
ity matrix elements for the components of the current 
(plane- wave values are shown following the arrows), 



h° = A 2X .cos- 



C 2Xi sin—t 



K = B 2X% sin^6 e + D 2Xi cos^i 



cos-V e , 
■ sin^6 e , 



h\ = 2iX i ( y A 2 \ i sin^6 e - C 2 \ z cos^6^j — > 2iXisin^6 e 



B 2Xi cos-6 



D 2 \ i sin—i 



cos-U e , 



(59) 



where the spin-orbit parts of the eikonal phases enter the 
helicity matrix elements in the following four combina- 
tions 

A 2Xi = cos^ f ~ ] * cos^ +) - sin^ +) sin^f^e 2 ^*'^ , 
B 2X% = cosj ( f ) * cos^ +) + sin^ +) sin^* e 2iX ^~' t >f\ 
C 2X% = cos*i { ~ ] * sin*/\ +] > e 2lX ^ + cos^ +) sin^f^e' 2 ^ , 
D 2X% = cos^f^sin^e 2 ^ - cos^sm^'e -2 '*'*'. 



general form expressed in terms of the complex phase 
7W = 7 < - ± - ) =p and based on the zero-mass limit for 
the lower components of the Dirac wave function. One 
may work at various orders of approximation by using the 
results of Section III for the eikonal phases. The leading 
order contributions are from 7p\ which is real. 

Some insight into the helicity matrix elements can be 
obtained by evaluating them for forward scattering and 
backward scattering. For forward scattering, the impact 
parameters and azimuthal angles for initial- and final- 
state eikonal phases are equal, i.e., bf = bi and 6f = 
(f>i. Evaluating the expression of Eq. (60) one finds the 
following simpler forms, 

A 2Xi =cos^ f - ) *+^ +) ), 

i?2A l = COS (7 4 (+) -7r ) *), 

C 2Xi =±sm(7}- ) *+7| +) ), 



D 2Xi =± S m( 7 j~ ) *-7| +) ) 



(61) 



where the upper sign applies for <j>i = 6f — and the 
lower one for 6i — 6/ — it. The helicity matrix elements 
for the same values of 6i are 



cos 



(^eT(7i +) +7^)), 



K - 8in(±0 e ± (7 4 (+) - l { f ] * 

/ i ^(2iA i )sin(^ e T(7! +) l7}" ), )) 1 

K ^cosQ-e e ±{^ -^ f - ) *)). (62) 

Note that the helicity matrix elements are similar to the 
plane-wave matrix elements except for a shift of the elec- 
tron scattering angle. The shift depends upon the az- 
imuthal angle Si and as indicated by the ± signs, the 
shifts at <pi = it are opposite to those at <pi = 0. Can- 
cellations are expected in the integration over 6i. These 
shifts that occur due to the spin-orbit interaction affect 
the longitudinal and transverse parts of the current in 
somewhat different ways. They may provide interesting 
insight into the accuracy with which one may make the 
L/T separation in the presence of Coulomb corrections. 
However, the required numerical evaluation is beyond the 
scope of this paper. 



QUASI-ELASTIC ELECTRON SCATTERING 
BY NUCLEI 



(60) 



Spin-orbit phases are of order 1/fc so the corrections 
that they produce are less important at higher electron 
energy. The current matrix element given above is a 



In this section we consider quasi-elastic scattering of 
electrons by nuclei but only taking into account the 
longitudinal current and the Coulomb corrections that 
arise from spin-independent terms in the eikonal expan- 
sion. The relevant matrix element involves one-photon 
exchange between the electron and a nucleon in the nu- 
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cleus and the cross section takes a well-known form, 



(63) 



where p is the momentum of the knocked-out nucleon and 
E p = \J M 2 + p 2 is its energy. The bar denotes an aver- 
age over initial helicities and a sum over final helicities. 
The quasi-elastic matrix element is 

M = S XfXi J d*r J ^e^-O-V**) x 

^(q,p)- (64) 



^(-) (r)/r ^ (r) ^ (r) ,___ 



?D{+), 



In the plane- wave impulse approximation (PWIA), 
Coulomb distortion of the electron waves is neglected so 
the integration over r produces J( 3 )(q — Q). The matrix 
element simplifies to 



M PWIA = g 



XfXi 



(65) 



where hp WIA denotes the helicity factors shown in 
Eq. (59) after the arrows. The PWIA cross section may 
be expressed in terms of longitudinal and transverse re- 
sponse functions, Rl and Rt, as follows 



da _ J Q 4 p , Q 2 1 p 



where 



and 



CMott = 



4a 2 Ejcos 2 ?f 

Q 3 



2Q 2 



(66) 



(67) 



(68) 



With Coulomb corrections included, the longitudinal 
matrix element of interest must take a gauge invariant 
form. This requires that the electron current must be 
conserved in the sense that 

J d 3 ^* (r) (u;j e - q ■ j e ) e-^^ (r) = 0, (69) 

and that the nuclear current should separately be con- 
served, 

J d 3 r¥->(r)e^ (cjJ% - q - J N ^(r) = 0, (70) 

where q is the photon three momentum. With Coulomb 
distorted waves, the photon momentum q differs from 
the electron's momentum transfer Q = k; — k/ and the 
longitudinal current is defined with respect to the direc- 
tion of the photon that is exchanged, not with respect to 



the difference of asymptotic electron momenta. Current 
conservation follows because the current obeys a Ward 
identity similar in form to Eq. (31). 

Owing to current conservation, the longitudinal cur- 
rent matrix element can be simplified as follows, 



M-(9-j e )(«-Jjv) = ieV 



jO i 



1 - 



(71) 



which holds either for the Klein-Gordon or Dirac case. 

When the Dirac equation is used for the electron's 
Coulomb distorted waves, the matrix element due to the 
longitudinal current is 



Mi 



(2n) 2 



/f ( - ) (r)/f (+) (r)/ l °(r)^Jj ^ ( q,p)(72) 

where Eq. (71) has been used to include the components 
of j e and Jjv that are parallel to q. The longitudinal re- 
sponse function is obtained by dividing the cross section 
integrated over the angles of the knocked-out nucleon by 
the Mott cross section, 



Rl = 



Q 4 



VMottQ 4 



J dn p -^\M L \ 2 E* P E pi (73) 



where Ml is the longitudinal amplitude of Eq. (72). The 
full calculation thus involves a six-dimensional integra- 
tion in order to obtain the amplitude Ml- Two more 
integrations over the angles of the knocked-out nucleon 
are required in order to obtain the response function. 
Results based on the eight-dimensional integration are 
called "full calculations" in the following sections. 

When the Klein-Gordon equation is used for the elec- 
tron's Coulomb distorted waves, the time component 
of the current is j° = [E, + E f - 2V{r)]/ y /AE l E f 
when Coulomb effects are included and j° = (E^ + 
Ef)/ y/AEiEf in the PWIA. The matrix element in the 
Klein-Gordon case is 



fr^' w* 1 ~l:' v/, [ -i K(q,p), (74) 



(2tt)' 

t KG{+)^ [Ej+E f -2V(r)} ( 1_ 

q 2 



V4£~E 



/ 



An analog of the Mott cross section based on the longi- 
tudinal current is 



rKG _a 2 E s {E i+ E } f 



EM A 



and the response function is 



Rl = 



1 



T KG 



|M£ G | EfpEp. 



(75) 



(76) 
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The longitudinal rcspsonse function Rl is calculated 
using a very simple model of the nuclear current as fol- 
lows, 



^(q,p) 



Pi+Pf 



y/4E p (E p - lu) 



U(q-P), (77) 



where ^»(k) is a gaussian wave function for a bound nu- 
cleon, 



$(k) = (27r/3 2 ) 3 / 4 e- 



(78) 



VI. EFFECTIVE MOMENTUM 
APPROXIMATION REVISITED 

As shown by Rosenfelder [25] and Traini [26], there 
are significant cancellations in the Coulomb corrections 
when response functions are evaluated in an effective- 
momentum approximation (emo). This approximation 
usually is based on expanding the eikonal phase in a Tay- 
lor's series about r = and keeping the first two terms 
as follows, 



X (r) «x(0) + r- V X (0 + - 



(84) 



normalized such that / d 3 fc|V>(k)| 2 /(27r) 3 = 1. This sim- 
ple model is used because the Coulomb corrections should 
depend mainly on the electron wave functions. Calcula- 
tions are based on the value (3 = 2 fermi. 

The nuclear current given above is based upon initial 
and final momenta, 



P>} = (E P7 p), 
Pi =(E p -u>, p q) , 



(79) 



where lu and q arc the photon's energy and momentum. 
Because of energy conservation, E p = M + lo — B, where 
B rts .008GeV is a typical binding energy of a nucleon. 
Gauge-invariance must hold so Eq. (70) is used to elimi- 
nate the component of the nuclear current that is parallel 
to the photon's momentum. 

In the plane-wave impulse approximation, the longitu- 
dinal response function is 



pE p 



Q 2 M L 



(2tt) 5 j - cos-_ 
where the longitudinal amplitude is 

^pwia^n (Q> P) 



(80) 



M PWIA 



= 5\ f > 



Q 2 



(81) 



Using the current and wave function described above 
leads to 



R 



PWIA _ P ( 2E P - 

(2tt)3 4(E p - cu) 



dO p |V(Q- P )| 2 



(82) 



The angular integrations are straightforward, yielding 

1 (2E p -lu) 2 13 



R 



PWIA 



(w,Q) 



2tt A{E p -uj) |Q| 



e -/3 2 (IQ|-p) 2 /2 _ e -/3 2 (|Q|+p) 2 /2 



(83) 



Here the PWIA response function is normalized so that 
at fixed Q, J cIluRl(Q,lu) w 1. 



The focusing factors are approximated by their values at 
r = and the helicity matrix elements are approximated 
by the plane-wave values. Integration over r then gives 
<j( 3 )(q _ Q e ^j), so the longitudinal amplitude simplifies 
to the PWIA form as follows, 



MT° = ^S XfXi h° PWIA J° N (Q eff ,p) 



ff ( -\o).fr +> w 



Q 2 



eff 



(85) 

The effective momentum involves the gradient of the 

eikonal phase shift \ — X/ ^ +Xi + ^ at the origin. Because 

of cylindrical symmetry of x\ + ^ about the direction fcj, 

Vx| + ' ) at the origin is nonzero only along the direction fcj, 

and similarly VxyT^ at the origin is nonzero only along 

the direction kf. With Vi = Vf w 1, we find the same 
result as Traini, 



Qeff = h hi — 5k 



k f k f 



5k 



(86) 



where 5k = V(0). It is correct up to first order in the 
eikonal expansion because the contribution from the gra- 
dient of eikonal correction xi vanishes at the origin. 

The analyses of Rosenfelder [25] and Traini [26] are 
based on the approximate focusing factors, //^ + ' w 

1 - V(0)/Ei and ff { ~ ] w 1 - V(0)/E f . Coulomb ef- 
fects in the focusing factors and the effective photon 
propagator cancel if one considers the photon propaga- 
tor of the transverse amplitude, which is l/[Q 2 ^j — to 2 ] — 
l/[4[k - V{0)][k f - V(O)]sm 2 ±0 e ], i.e., 



,£>(- 

h 



(0)/f (+) (0) 



Q 



eff 



- LU^ 



1 



(87) 



which is the same as in the plane- wave case, Eq. (65). 
However, the requirements of gauge invariance that have 
been incorporated into the longitudinal matrix element 
of Eq. (72) show that the effective photon propagator is 
l/Q 2 y-y. In that case the cancellations still are significant 
but not perfect. We find 



/f (_) (0)/, 



D(- 



'(0) 



Q 2 



eff 



Q 2 + A' 



(88) 
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where A represents a Coulomb correction to the plane- 
wave result, 



2 V(0) 



[hi + k f - V(0)} 
[h-V(0)}[k f -V(0)} 



(89) 



As an example, for 500 MeV electrons scattering from a 
25 MeV Coulomb potential with # e =60°, A/Q 2 w .01 at 
uj = .125, which is close to the quasi-elastic peak. 

Numerical calculations based on partial-wave expan- 
sions of Dirac- Coulomb waves also indicate that the 
Coulomb corrections do not cancel to the extent that 
Eq. (87) would suggest. [43] For the e + and e— response 
functions the ema analysis suggests that the same re- 
sponse function should be obtained if the energy is shifted 
such that Ei(e+) + V(0) = Ei{e~) - V(0), where V(0) is 
the electron-nucleus potential at r = 0. This shift results 
in the same Q e // for e+ and e— scattering. Using the 
gauge-invariant response function, the correction A has 
opposite sign for e + scattering than for scattering. It 
can produce a 4% difference in the e+ and e~ longitudi- 
nal response functions, whereas there would be no differ- 
ence based on Eq. (87) . Numerical calculations based on 
partial- wave expansions of Dirac-Coulomb waves indicate 
that the Coulomb corrections do not cancel to the extent 
that Eq. (87) would suggest. In their analysis based on 
the full DWBA, Kim et al. [44] have obtained results for 
the sum of longitudinal and transverse responses, which 
differ by about 15% for 420 MeV e+ and 383 MeV e-. 



VII. NUMERICAL CALCULATIONS FOR 
QUASI-ELASTIC SCATTERING 

Calculations of the longitudinal response function are 
performed for three cases: PWIA, ema and the full cal- 
culation. We also use distorted waves based on the 
Dirac equation and the Klein-Gordan equation. The 
calculations are based on a charge Z = 38 and ra- 
dius R = 2 fermi in the Coulomb potential of Eq. (48). 
Eikonal phases are evaluated through second order, i.e., 
X = Xo + Xi + X2 and u> = lj\ + LU2- However, the expan- 
sion converges rapidly for the parameters and energies 
used and results based on xo + Xi an d oji differ by about 
0.3% at the quasi-elastic peak. 

In order to check the accuracy of the full calculation of 
Rl, which involves an eight-dimensional integration, we 
have performed calculations with Z = and compared 
with the analytical PWIA result. Results are accurate 
to between 0.5% to 2% with the integration points that 
have been used. 

Figure 2 and Table I show the longitudinal response 
function for 500 MeV electrons with scattering angle 
e = 60° based on the use of the Dirac current and dis- 
torted waves. In Fig. 2 the full calculation is shown by 
circles, the ema calculation based on 5k = V(0) is shown 
by the solid line and the PWIA calculation is shown by 
the dotted line. Figure 3 shows the response function 
for the same kinematics using the Klein- Gordon current 




0.2 

(GeV) 



FIG. 2: Longitudinal response function versus the electron's 
energy loss, ui, calculated using Dirac distorted waves for e~ 
scattering at E = 500MeV and 9 e = 60°. Dotted line shows 
PWIA, solid line shows ema and circles show full calculations 
based on Eq. (72). 



TABLE I: Numerical results for R L 
scattered by 60°. 



for 500 MeV electrons 



GeV 


PWIA 


ema 


Full 


0.040 


0.49 


0.23 


0.24 


0.050 


1.22 


0.65 


0.73 


0.060 


2.41 


1.42 


1.62 


0.070 


3.98 


2.59 


2.92 


0.080 


5.73 


4.07 


4.55 


0.090 


7.35 


5.68 


6.20 


0.100 


8.56 


7.16 


7.63 


0.110 


9.18 


8.28 


8.65 


0.120 


9.16 


8.87 


9.12 


0.130 


8.57 


8.90 


8.84 


0.140 


7.59 


8.42 


8.24 


0.150 


6.39 


7.55 


7.27 


0.160 


5.14 


6.45 


6.12 


0.170 


3.98 


5.29 


4.95 


0.180 


2.96 


4.17 


3.85 


0.190 


2.14 


3.17 


2.90 


0.200 


1.50 


2.34 


2.12 


0.210 


1.02 


1.68 


1.51 


0.220 


0.68 


1.17 


1.06 



and distorted waves. The results are quite close to those 
based on the Dirac equation, which is expected because of 
the similarity of eikonal wave functions for the two cases 
when the spin-dependent phases are omitted. As noted 
earlier, the Klein-Gordon and Dirac results involve simi- 
lar overall focusing factors when consistent currents and 
wave functions are used. We also show by the dashed 
line in Fig. 3 a calculation following the prescription 
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of Ref [35] in which the Klein- Gordon current is used 
with the Dirac focusing factors. This produces a signif- 
icantly different result because of the extraneous factor 
1 — V(0)/E that is included. The reason why the extra 
factor produces a large change is because the other fo- 
cusing factors largely cancel out of the matrix element as 
in Eqs. (87) or (88).' 
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FIG. 3: Longitudinal response function versus the electron's 
energy loss, u), calculated using Klein-Gordon distorted waves 
for e" scattering at E = 500MeV and 6 e = 60°. Dotted 
line shows PWIA, solid line shows ema, dashed line shows 
the prescription of Ref. [35] and circles show full calculations 
based on Eq. (74). 

For the remainder of this section we use only the Dirac 
current and Dirac wave functions. Figure 4 shows the 
longitudinal response function for e + scattering at 540 
MeV. In general the ema is seen in Figures 2 and 4 to 
produce a significant shift of Rl away from the PWIA 
result and towards the full calculation of R L . There also 
is good agreement between the response functions for e— 
and e+ scattering at the energies that make Q e // close 
to the same for both. 

Upon closer inspection, we find that the ema result 
is not precise. Figures 5 and 6 show the ratios of the 
ema response functions to the full ones for e~ and e + 
scattering, respectively. Deviations of 5-10% occur at 
values of u> that are away from the quasi-elastic peak. 
In the literature one finds a variety of suggestions for 
improving the ema, such as using 5k = V(R) rather than 
5k = V(0), where R is the mean nuclear radius, or using 
5k = V, where V is the average potential within a sphere 
of radius R. These prescriptions improve the ema at 
some values of w but not all values. 

A precise form of the effective-momentum approxima- 
tion would be useful for removing Coulomb corrections 
from experimental data in a straightforward manner. It 
would allow a determination of the PWIA response func- 
tion. In order to have a precise result, one should de- 
termine appropriate values of the momentum-shift func- 
tion 5k(ki,u),9 e ) from which the appropriate Q e // may 
be calculated as in Eq. (86). In order to determine this 
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FIG. 4: Longitudinal response function versus the positron's 
energy loss, w, for e + scattering at E — 540MeV and 9 e = 
60°. Dotted line shows PWIA, solid line shows ema and 
circles show full calculations based on Eq. (72). 
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FIG. 5: Ratio R L ma / Rl as a function of energy loss u> for 500 
MeV e~ scattering at angle 9 e - = 60°. 



function, we have fit the full response function as follows, 
R L (Q,w) = A Rl WIA (Q eff ,w) (90) 



at fixed values of fc, and 8 e by varying 5k and the nor- 
malization constant A. The role of the normalization 
constant A is to ensure that 5k(u>) is a smooth function. 
Without this parameter, 5k can take anomalous values 
near the quasi-elastic peak. Figures 7 and 8 show the val- 
ues of 5k that are found to yield Eq. (90) for scattering of 
e" (with A e - = 0.987) and e+ (with A e + = 0.973). The 
sign of 5k is in accord with the potential V(0), which is 
-27 MeV for electrons and +27 MeV for positrons. How- 
ever, a constant value of 5k does not suffice. The required 
values vary with u> at a fixed value of 8 e = 60°, within a 
range 0.7|V(0)| < \5k\ < 1.5|V(0)|. Note that the shift 
5k required for the e + response function is not simply 
obtained by reversing the sign of the shift required for 
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FIG. 6: Ratio R e L ma /R L as a function of energy loss w for 540 
MeV e + scattering at angle 8 e = 60°. 
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FIG. 8: Shift 5fc used in Q e // in order to achieve the equality 
of Eq. (90), as a function of positron's energy loss, lu, for 540 
MeV e + scattering at 9 e + = 60°. The normalization constant 
is A = 0.973. 
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FIG. 7: Shift 8k used in Q e // in order to achieve the equality 
of Eq. (90), as a function of electron's energy loss, u>, for 500 
MeV e~ scattering at 6 e - = 60°. The normalization constant 
is A = 0.987. 



the e~ response function. Nevertheless, the closeness of 
A to unity and 5k to V(0) suggests that an ema analy- 
sis using the fit values of 5k would be well-motivated on 
physical grounds. 

If the momentum-shift function and normalization con- 
stant A are determined theoretically for a given nucleus 
based on a sophisticated model of the nuclear current, 
and Rl values are available based upon experimental 
data, our results suggest that one may use Eq. (90) 
with the experimentally determined Rl in order to ex- 
tract information about the undistorted response func- 
tion, R^ WIA . Of course, the accuracy would depend 
upon the accuracy of the L/T separation. 



VIII. SUMMARY 

In this paper, we develop the eikonal expansion for 
relativistic wave functions based on the Klein-Gordon 
equation and the Dirac equation, each with a Coulom- 



bic potential. The purpose is to obtain some insight into 
the Coulomb corrections in quasi-elastic electron scatter- 
ing without a significant loss of accuracy. The eikonal 
expansion is carried out to obtain corrections up to or- 
der 1/fc 2 to the eikonal approximation. We show that 
focusing factors are obtained in a systematic manner by 
use of the eikonal expansion. Although focusing factors 
take somewhat different forms for the Klein-Gordon and 
Dirac wave functions, equivalent results are obtained for 
the current matrix elements for the two cases because of 
the Coulomb correction to the Klein-Gordon current. 

Based on a simple form of the Coulomb potential, an- 
alytical results are given for the eikonal phases and thus 
for the eikonal wave functions. For scattering of electrons 
with energies of a few hundred MeV or more, the eikonal 
expansion converges rapidly. 

Coulomb corrections in quasi-elastic electron scatter- 
ing are considered using the analytical wave functions 
based on the eikonal expansion. For the longitudinal re- 
sponse function, we show that the approximate evalua- 
tion of the matrix element using the effective-momentum 
approximation is modified by the requirements of gauge 
invariancc. That modification causes a small but signif- 
icant difference between the e~ and e + response func- 
tions. Moreover, the ema is not sufficiently accurate to 
allow a precise analysis of data because the effective mo- 
menta ki — 5k and kf — 5k are not precise when 5k is 
taken to be independent of u>, as it is in the usual form 
of ema. Using a simple model of the nuclear current, 
we find that use of a function 5k(ui) can yield a precise 
form of the ema. The analysis should be repeated for 
more sophisticated models of the nuclear current with 
the goal of determining the function 5k(uj) for different 
models. Based on our results, we suggest that if one is 
able to extract Rl(w, E,8 e ) from experimental data at 
fixed electron beam energy, E, and fixed electron scat- 
tering angle, 9 e , then it may be equated to a constant 
A k 1 times the PWIA response function evaluated at 
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an effective momentum transfer. 

Whether one may extract the longitudinal response 
with reasonable accuracy has not been resolved in this 
work. We have addressed the transverse current ma- 
trix elements and have shown how the spin-dependent 
Coulomb corrections modify the plane- wave helicity ma- 
trix elements. There is a potentially interesting effect 
in the manner that Coulomb corrections enter as a shift 
of the electron scattering angle in helicity matrix ele- 
ments. However, these spin-dependent Coulomb correc- 
tions have not been included in our calculations. Eval- 
uating their effects may help to provide insight into the 
accuracy with which longitudinal and transverse response 
functions can be separated. Because we have only calcu- 
lated the response functions based on the longitudinal- 



current matrix elements, our results cannot be compared 
with experimental data, or with a full DWBA analysis, 
because both include an inseparable sum of longitudinal 
and transverse responses. In particular, we note that the 
effects of Coulomb corrections in the helicity matrix ele- 
ments of Eq. (62) would go in different directions for e+ 
and e— scattering. 
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APPENDIX A: ANALYTICAL PHASES 

Analytical results are given for the cikonal phase shifts 
based on the potential of Eq. (48). Terms of order m 2 /E 2 
are omitted. In order to give the analytical results in 
a reasonably short form, we define the following quan- 
tities: r\ — Za/v, u — \Jr 2 + R 2 , w = \/b 2 + R 2 and 
/3 = tan~ 1 (w/z), where [3 varies from to 7r as z varies 
from — oo to oo. In Xa^\ a term r/ln(2A) is omitted, 



where the eikonal integral is calculated over the range 
{— A, z}. This contribution provides an irrelevant phase 
factor that occurs for potentials that behave as 1/r as 
r — > oo. 

The analytical results are 



Xo +) = nin 



z + u 



Xi +) 



V 2 b 2 
kw 4 



z + u 



(Al) 



(A2) 
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APPENDIX B: EIKONAL EXPANSION FOR 
DIRAC WAVE FUNCTION 

For z — > — oo, the incoming wave should reduce to a 
plane wave, which is realized by the boundary conditions 



Eq. (34) leads to 



as z — > — oo. Inserting Eq. (35) into 
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The quantity V 2 Xo + ^ that occurs in x^ f° r tne Dirac 
case is given by 



v 2 x^ +) = v 4 

u 6 



z + 2u Au 3 



(z + u) 2 



ur 



(A9) 



That term should be omitted for the Klein-Gordon case. 

Corresponding results for the incoming wave boundary 
condition are obtained from Eq. (24), which holds for any 
of the quantities in Eqs. (Al) to (A9). 



Momentum operators in this expression act on all quan- 
tities to their right. Performing the indicated differen- 
tiations of the factor (1 — V/E 2 ) 1 ^ 2 and multiplying by 
El ,2 {E 2 - V) 1 ' 2 /{2k) leads to 

( dx (+) Vc (Vx (+) ) 2 , V 2 X (+) «eV^ +) 
\ dz v 2k 2k 2k 



. d .Vx (+) W 
+l d^ + l ~k~ ' " 2k(E 2 - V) 



• V 



+ 2k(E 2 -v) v + 2k ) 



(B2) 



where v = k/E 1 and we have defined central and 
spin-orbit potentials as follows, 



V 2 (r) V 2 V(r) 



3(W(r)) 5 



and 



2E AE(E 2 - V(r)) 8E(E 2 - V(r)) 2 ' 
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db E 2 -V dz 



db 



The spin-orbit term in Eq. (B2) involving potential vi + ' 
gives rise to the spin-dependent eikonal phase 

Some care is required in evaluating the derivatives of 
the spin-dependent phase factor, for example, Vcr e = 
-{ub/b)e and V 2 cr e = -cr e /b 2 . We find 



d X (+) V c (Vx (+) ) 2 , .V 2 x (+) a e V t 

\ 1 
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dz v 



2k 
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2k 
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where an overall factor e l,Te7<+> is omitted, and 



2fc(£ 2 - V) 



(B6) 



The quantity is found to be of order 1/fc 3 or smaller 
so will be dropped in the following. 

The part of Eq. (B5) that involves a e and the part that 
does not involve a e must vanish separately. Using trace 
techniques to project out these parts produces equations 
for the phases and 7^ as follows, 



dx {+) 

dz 
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One may see from the second of these equations that 7 ( + ) 
is of order 1/fc or smaller. Thus, the last two terms in 
Eq. (B7), the last term in Eq. (B8) and X 3 are of order 
1/fc 3 or smaller. We omit these terms in the following 
with the objective of obtaining results correct to order 
1/fc 2 . We also simplify sin(2 7 ( + )) w 2 7 (+) for the same 
reason. 

Integration in accord with the outgoing-wave boundary 
conditions produces the basic equations 



;(+)(r) 



dz 



,V C {V) 



dz 



,(V'x (+) (r')) 2 



2fc 



J — ( 



+i / dz 

-00 



,V 2 X (+) (r') 



2fc 



(B9) 



and the complex phases \( \r) and 7 ' )(r) obey 
^(-) (r) = _ 
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where vj ' is 



1 dV(r) 9x (_) (r) 



£ 2 - V dz 



db 



(B14) 



The eikonal expansion is the iterative solution appro- 
priate to large k of Eqs. (B9) and (BIO) or Eqs. (B12) 
and (B13). The expansion for x initially takes the 
form 



x (+) =x^ +) +xi +) +x 2 +) + -' 



(B15) 



where each term in the series is smaller than the previous 
one. The "barred" phases are in general complex with 
the exception that Xo + ^ is rea l when the potential is real. 

Lowest order terms may be obtained from Eqs. (B9) 
and (B12) by keeping just the leading terms on the right 
side, 



(B16) 
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where r' = (b, z'). 

With incoming- wave boundary conditions, the wave 
function must become a plane wave as z ^ +00 and 
is written as 



(l-_ 



V 



1/2 



e ikz e -ix { ' e -iff e 7 



v(-) 



• (BH) 



It should be noted that the same symbol Xo was used 
in the Klein- Gordon case but the meaning is different 
here because of the extra terms in V c (r), Eq. (B3). Us- 
ing the same iterative procedure that was used to obtain 
Eqs. (10) and (11) leads to 

x! +) (r) = ~ fjz'(W'xi +) (r')) 2 

+ Yk fj^' 2 xt\*'l (B17) 

X W(r) = -i- fjz'[2V'x[ +) (v')-Vx[ + \v') 
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For the spin-dependent phase, 7^ + \ the iterative ex- 
pansion is 

7 (+) =# ) +# ) +- (B19) 

where 

7r +) (r) = -±- J Z J z ' V (+)(h,z'), (B20) 



Complex eikonal phases and are decomposed 
into real and imaginary parts ( the "unbarred" phases) 
as follows 

X (+) = x (+) + ^ (+) 

7<+) = 7 (+)+ij(+). (B22) 

Because the central and spin-dependent potentials have 
terms that involve various powers of 1/fc, the "barred" 
phases that are obtained from the iterative expansion 
above are not ordered systematically. Therefore a second 
expansion is performed to obtain terms that are propor- 
tional to powers of 1/fc. This leads to Eqs. (36) to (37) 
of the text. 



